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ABSTRACT 

Utilization  of  composites  in  critical  design  applications  requires  an 
extensive  engineering  experience  data  base  which  is  generally  lacking, 
especially  for  rapidly  developing  constituent  fibers.  As  a  supplement,  an 
accurate  reliability  theory  can  be  applied  in  design.  This  investigation  is  a 
part  of  a  research  effort  to  develop  a  probabilistic  model  of  composite 
reliability  capable  of  using  data  produced  in  small  laboratory  test  samples  to 
predict  the  behavior  of  large  structures  with  respect  to  their  actual 
dimensions.  This  work  included  testing  of  composite  strength  which  was 
then  used  in  exploring  the  methodology  of  predicting  composite  reliability 
from  the  parent  single  filament  fiber  strength  statistics.   This  required 
testing  of  a  coordinated  set  of  test  samples  which  consisted  of  a  composite 
and  its  parent  fibers.   Previously  collected  fiber  strength  statistics  from  two 
different  production  spools  were  used  in  conjunction  with  the  current  effort. 
This  investigation  established  that,  for  a  well  made  composite,  the  Local 
Load  Sharing  Model  of  reliability  prediction  exhibited  outstanding 
correlation  with  experimental  data  and  was  sufficiently  sensitive  to  predict 
deficient  composite  strength  due  to  a  specific  fiber  spool  with  an  abnormally 
weak  lower  tail.  In  addition,  it  provided  an  upper  bound  on  the  composite 
reliability.  This  investigation  is  unique  in  that  it  used  a  coordinated  set  of 
data  with  an  unambiguous  genesis  of  parent  fiber  and  subsequent  composite. 
The  findings  of  this  investigation  are  also  definitive  in  that  six  orders  of 
extrapolation  of  size  in  reliability  prediction  has  been  verified. 
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I.     INTRODUCTION 

As  more  and  more  refinements  have  been  developed  on  composite 
material  during  the  past  several  years,  the  use  of  these  materials  in 
structural  applications  has  become  commonplace,  if  not  the  material  of 
choice,  in  industry  today.  The  reasons  for  this  are  well-known,  the  most 
common  being  the  excellent  strength-to-  weight  and  stiffness -to- weight 
ratios  associated  with  many  composite  materials.  In  addition,  through 
micro-mechanical  analysis,  the  engineer  can  design  the  material  to  optimize 
the  physical  properties  (such  as  strength,  thermal  or  hygro  properties)  for 
the  functional  requirements  of  the  structure.  As  in  any  new  developmental 
process,  improvements  to  the  product  are  always  sought  and  the  recent  past 
has  seen  rapid  growth  in  fiber  improvements  with  respect  to  fiber  strength 
characteristics.    These  new  and  improved  fibers  are  then  being  used  in  the 
production  of  todays  composite  material.  Associated  with  the  rapid  growth 
in  fiber  improvements  is  the  general  lack  of  an  experience  data  base 
regarding  the  reliability  and  maintainability  of  the  structure  throughout  its 
operational  life. 

The  Advanced  Composites  Laboratory  at  the  Naval  Postgraduate 
School,  Monterey,  CA,  is  addressing  these  questions  with  a  systematic 
approach  to  composite  reliability  research.  Unlike  traditional  design  and 
manufacturing  processes  utilizing  materials  with  well  known  characteristics 
(such  as  aluminum)  where  very  large  design  "data  bases"  are  known  which 
often  permits   "design  by  experience",  the  incorporation  of  composite 


materials  in  design  and  manufacturing  is  still  in  its  infancy,  thus  an 
equivalent  "data  base"  does  not  exist.  Design  of  a  composite  structure  must 
be  based  on  a  sound  reliability  model  of  the  salient  physical  failure  process. 
Under  tension,  the  failure  process  of  a  fiber  composite  is  sequential. 
Internal  failures  start  with  breaking  of  weak  fibers  within  the  composite  at 
relatively  low  applied  loads.  The  original  stresses  carried  by  the  broken 
fibers  are  transferred  by  the  matrix  binder  to  be  shared  by  the  neighboring 
unbroken  fibers.  This  load  sharing  mechanism  provides  local  redundancy 
thereby  delaying  global  catastrophic  failure.  With  the  addition  of  external 
loads,  more  fiber  failure  sites  are  created  leading  to  clustering  of  the 
failure  sites  which  leads  to  ultimate  catastrophic  failure.  This  failure 
process  is  modeled  in  the  Local  Load  Sharing  model.  Also,  questions 
dealing  with  size  effect  (strength  is  inversely  proportional  to  physical 
dimensions)  of  the  composite  structure  must  be  addressed  during  the  design 
process.   This  investigation  is  a  part  of  the  research  effort  to  develop  a 
probabilistic  model  of  composite  reliability  which  is  capable  of  predicting 
actual  large  dimension  structures  from  data  produced  in  small  dimension 
laboratory  samples. 

Composite  reliability  is  a  function  of  both  the  strength  and  life  of  the 
constituent  fiber  and  matrix  material.  To  establish  a  probabilistic  model  of 
a  composite's  reliability,  an  adequate  number  of  experiments  must  be 
conducted  which  deal  with  the  strength  and  life  of  not  only  the  composite 
but  also  of  its  parent  constituent  fibers.  The  objective  of  this  study  was  to 
gather  composite  strength  reliability  test  data  to  be  used  in  exploring  the 
viability  of  predicting  composite  reliability  using  a  Local  Load  Sharing 


(LLS)  model  produced  from  previously  obtained  single  fiber  strength 
statistics.    This  was  accomplished  by  testing  composite  strands  to  failure 
and  comparing  the  data  with  the  LLS  reliability  prediction. 


II.      BACKGROUND 


A.  FIBER  REINFORCED  COMPOSITE  STRAND  FAILURE 
MODEL 

Graphite  composite  materials  usually  consist  of  high-strength,  high- 
stiffness  fiber  encased  in  a  ductile  matrix.  The  strength  characteristics  of 
the  composite  are  dependent  on  the  modulus  of  the  fiber,  the  modulus  of 
the  matrix  and  the  effectiveness  of  the  bond  between  the  matrix  and  the 
fibers  at  their  interface.  The  failure  model  used  in  this  study  was  first 
introduced  by  B.  Walter  Rosen  [Ref.  1].  The  model  addresses  the  failure 
of  a  composite,  consisting  of  a  matrix  stiffened  by  uniaxial  oriented  fibers 
when  subjected  to  a  uniaxial  tensile  load  parallel  to  the  fiber  direction.  It 
assumes  that  the  fibers  have  a  statistical  distribution  of  flaws  or 
imperfections  which  results  in  individual  fiber  failures  in  the  composite  at 
various  stress  levels.  The  composite  itself  will  fail  when  the  remaining 
unbroken  fibers,  at  the  cross-section  with  the  most  fiber  failures  (weakest 
cross-section),  are  unable  to  carry  the  applied  load.  Therefore,  in  this 
model,  the  composite  failure  is  dominated  by  the  tensile  fracture  of  the 
fibers.  In  analysis  of  the  statistical  distribution  of  flaws  in  an  individual 
fiber,  it  has  been  shown  [Ref.  2  and  3]  that  the  fiber  segments  can  be 
modeled  as  a  series  of  links  with  the  fiber  as  a  chain.  The  chain  (fiber) 
will  fail  when  the  applied  load  is  large  enough  to  cause  failure  of  the 
weakest  link.  The  statistical  strength  distribution  of  the  links  can  therefore 


be  approximated  by  extreme  value  distributions  which  may  be  expressed  as 
a  two  or  more  parameter  Weibull  distribution.   Such  distributions,  when 
represented  in  a  Weibull  Probability  of  failure  plot  appear  as  linear  or 
piecewise  linear  lines. 

If  one  were  to  combine  many  individual  fibers  into  a  bundle  and  a 
uniaxial  load  was  applied  to  it,  individual  fibers  would  begin  to  fail  as  the 
load  on  each  fiber  exceeded  that  of  its  weakest  link.  This  initial  failure 
load  is  much  less  than  that  of  the  bundle  failure  load.  Once  an  individual 
fiber  has  failed,  the  load  that  it  was  carrying  is  now  distributed  over  the 
remaining  fibers  within  the  bundle,  increasing  the  applied  load  on  each 
individual  fiber.   The  broken  fiber  no  longer  carries  any  of  the  load.   If 
the  combination  of  the  applied  and  the  assumed  load  on  the  individual 
fibers  remaining  is  less  than  its  weakest  link  strength,  then  the  fiber  will 
remain  intact.   If  it  is  greater,  then  the  fiber  will  fail  and  the  same  transfer 
of  loading  takes  place.  This  process  will  continue  as  the  applied  load  is 
increased  and  bundle  failure  will  occur  once  every  individual  fiber  has 
been  broken. 

When  the  fiber  bundle  is  encased  in  a  matrix,  the  failure  process  of  the 
newly  formed  strand  is  changed  from  that  of  the  bundle.  The  model  that 
describes  the  strand  failure  consists  of  parallel  fibers  in  an  otherwise 
homogeneous  matrix.  Its  fiber  statistical  strength  distribution  is  assumed  to 
be  the  same  as  before.  Figure  1,  Appendix  A  shows  the  failure  model 
presented  by  Rosen  [Ref.l].  It  should  be  noted  that  in  this  model,  the 
extensional  stresses  in  the  matrix  are  neglected  relative  to  those  in  the  fiber 
and  the  shear  strains  in  the  fiber  are  neglected  relative  to  those  in  the 


matrix.   This  is  a  reasonable  assumption  for  fibers  that  are  very  strong  and 
stiff  relative  to  the  matrix.  As  a  uniaxial  load  is  applied  to  the  strand,  the 
fibers  with  the  weakest  links  begin  to  fail  just  as  described  in  the  bundle, 
but  unlike  the  bundle,  the  matrix  provides  a  unique  load  sharing 
characteristic  that  enables  the  internal  broken  fiber  to  continue  carrying  an 
applied  load.  In  the  vicinity  of  an  internal  failure,  the  uniaxial  load  that 
was  carried  by  the  fiber  is  transmitted  by  shear  through  the  matrix  to 
adjacent  fibers.  These  adjacent  fibers  now  carry  a  load  increased  by  a  load 
concentration  factor  (Kr).  The  internal  failures  result  in  shear  stresses  that 
locally  may  attain  very  high  values.  At  the  fiber  end  (point  of  failure), 
stress  (o)  goes  to  zero  and  the  shear  (x)  goes  to  a  maximum  value. 
Therefore,  a  portion  of  the  fiber  near  the  fiber  end  is  ineffective  in 
resisting  the  applied  load.  This  distance  is  known  as  the  ineffective  length 
(5).  Rosen  has  shown  [Ref.  1]  that  the  ineffective  length  (normalized  to  the 
fiber  diameter)  of  this  model  to  be 
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where         6f    is  the  fiber  diameter 

vf    is  the  volume  fraction  of  fiber  in  the  composite 
Ef    is  the  modulus  of  the  fiber 
Gm  is  the  shear  modulus  of  the  matrix 

({>     is  the  fraction  of  the  undisturbed  stress  value  below  which 
the  fibers  shall  be  considered  ineffective 


As  the  applied  load  increases,  the  number  of  statistically  distributed  fiber 
failures  increase,  producing  an  accumulation  of  ineffective  fiber  lengths. 
When  a  sufficient  number  of  these  ineffective  fiber  lengths  combine  in  the 
vicinity  of  one  cross-sectional  area  of  the  strand,  it  results  in  a  weak 
surface.  This  enables  the  onset  of  matrix/fiber  interface  debonding  and/or 
crack  propagation  through  the  matrix,  ultimately  resulting  in  the  composite 
failure.   In  summary,  the  models  treat  fiber  failures  as  the  result  of  a 
statistically  distributed  flaw  (weakest  link).    Composite  strand  failures  are 
the  result  of  a  statistical  accumulation  of  the  fiber  failures  over  a  given 
cross-sectional  area.   Therefore,  fiber  strength  is  dependent  on  length;  that 
is,  longer  chains  (fibers)  have  a  higher  probability  of  having  a  weaker  link 
(flaw)  than  a  shorter  chain.  This  agrees  with  experimental  data  [Ref.  4] 
that  demonstrates  that  fiber  strength  is  a  monotonically  decreasing  function 
of  fiber  lenuth. 
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B.  LOCAL  LOAD  SHARING 

The  Local  Load  Sharing  model  used  in  this  study  was  developed  by 
D.  G.  Harlow  and  S.  L.  Phoenix  [Ref.  2  and  3]  and  is  shown  in  Figure  2, 
Appendix  A.  The  fiber-matrix  composite  is  viewed  as  a  planar  structure  of 
n  parallel  fibers  partitioned  into  a  series  of  m  sections,  called  bundles,  each 
with  n  fibers.  Conceptually,  each  bundle  (defined  as  the  metric)  can  be  any 
desired  length.  Harlow  and  Phoenix  designated  the  metric  length  to  be  the 
ineffective  length  8.  The  length  of  the  material  is  then  simply  the  metric 
length  times  m.  In  the  model,  the  bundles  are  considered  to  be  statistically 
independent  and  the  strength  of  the  composite  material  is  that  of  its  weakest 


bundle.  In  addition,  it  is  assumed  that  the  strength  of  the  fiber  elements  are 
statistically  distributed  and  can  be  modeled  by  the  Weibull  distribution 


F(L)  =  1  -  EXP  \  -  (%)   f         for  L>0 


where  a    is  the  shape  (slope)  parameter 

P    is  the  location  parameter 

L    is  the  independent  variable 
As  shown  in  Figure  2,  Appendix  A,  if  the  bundle  load  is  x  (per  fiber)  then 
the  adjacent  surviving  fiber  elements  carry  a  new  load  Krjc  where 

Kr  =  1  +r/2         r=  1,2,3 

and  r  is  defined  as  the  number  of  consecutive  failed  fiber  elements 
immediately  adjacent  to  the  surviving  element  (on  both  sides).  Kr  is  called 
the  load  concentration  factor.  The  probability  distribution  of  the  strength 
of  a  composite  with  Weibull  distributed  parent  fibers  can  be  computed 
numerically  using  a  recursive  relation  for  different  fiber  failure  cluster 
(Kr)  configuration.  Based  on  this  exact  calculation,  an  approximate 
representation  is  available  [Ref.  5]  and  is  summarized  as  follows: 


mn 


Hmn(L)=l-(l-W(L)r         forL>0 


where  W(L)  =  minF(k)(L)  for  L>0  and  k=  1,2, 3. 


and  F(k)(L)=l-EXP\-dk(1/n)     /  forL>0 


dk  =  dk(a)  =  2       (KjK^ K^) 

Kr  =  l+r/2       where  r  =  1,2,3 


with  the  variables  Hmn(L)  =  Probability  of  failure  of  a  composite  of  length 

m  times  the  metric  length 

W(L)  =  Probability  of  k  -  failure  in  the  bundle 

F(L)  =  Probability  of  failure  of  the  mn  fiber 
elements 

k    =  The  number  of  adjacent  fiber  breaks  in  a 
bundle 

dk  =  A  constant  that  is  a  function  of  a 

and  Kr,  r,  a,  (3  and  L  are  as  previously  defined.   Finally,  it  is  important  to 

note  that  for  a  poorly  made  composite  with  inferior  fiber/matrix  adhesion 

and  voids,  the  failure  process  will  approach  that  of  the  global  load  sharing 

of  a  bundle  rather  than  the  local  load  sharing  of  a  composite. 


9 


III.        METHOD  OF  TESTS 

A.  DESCRIPTION  OF  TEST  EQUIPMENT 
1.     Mechanical  System 

The  strand  tensile  strength  testing  was  conducted  utilizing  an 
tnstron  Universal  Testing  Instrument,  model  4206,  shown  in  Figure  3, 
Appendix  A.  It  is  comprised  of  two  major  systems:  a  crosshead  drive  and 
control  system,  which  applies  tensile  or  compressive  loading  to  a  specimen; 
and  a  highly  sensitive  load  weighing  system,  which  measures  the  load 
applied  to  a  specimen. 

An  Instron  2512  series,  1000  Kilogram  load  cell  which  was  used  to 
measure  applied  loads,  was  mounted  in  the  moving  crosshead  which  is 
operated  by  two  vertical  leadscrews  within  the  loading  frame. 

The  test  articles  (composite  strand  samples  physically  resembling  a 
pencil  lead)  were  held  by  two  Instron  Modular  Hydraulic  Wedge  Collet 
Grips  (series  2742),  one  mounted  on  the  load  cell  and  the  other  on  the  base 
of  the  machine.  One-half  inch  diameter  collets  were  used  in  the  wedge 
grips  for  mounting  the  specimen  into  the  grip.  The  specimens  were  held 
stationary  in  the  grips  by  hydraulic  pressure  obtained  by  a  Instron  Electric 
Pump  Pressurization  System  (model  A7154,  Rev  D,  #5215)  shown  in 
Figure  3,  Appendix  A.  The  unit  is  self-contained  with  gripping  pressure 
fully  adjustable  and  independent  of  system  pressure  or  fluctuations  such 
that  gripping  force  remains  constant  on  the  sample  under  test. 
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2.     Control  and  Data  Acquisition  System 

The  control  and  data  acquisition  of  the  testing  instrument  was 
automated  with  the  use  of  Instron  Series  IX  Materials  Testing  System.  This 
system  is  a  software  package  that  interfaces  the  test  instrument  with  a 
personal  computer.  The  testing  system  components  included  the  Series  IX 
software,  an  IBM/PC- AT  computer  (1.5  MB  of  RAM,  30  MB  hard  disk) 
with  an  EGA  monitor,  an  IEEE-488  (GPIB)  interface,  an  IBM  Proprinter 
II  parallel  printer  and  an  HP  7470A  plotter  for  graphic  functions. 

The  software,  which  is  menu  driven,  presents  the  operator  several 
options  including  computer  controlled  testing,  reanalyzation  of  data, 
calibration  of  the  test  instrument,  creating  and  or  modifying  the  test 
method  to  be  used  in  the  computer  controlled  testing,  and  plotting  of  the 
raw  data.   Separate  test  methods  were  created  for  each  gauge  length  of 
samples  tested  and  are  presented  in  Figures  4,  5,  6  and  7,  Appendix  A. 

B.  TEST  SPECIMENS 

1.     Sample  Composition 

The  test  specimens  were  produced  from  a  Hercules  Magnamite  high 
strength  graphite,  type  AS -4.  Both  the  fiber  samples  and  the  composite 
samples  were  fabricated  from  two  specific  production  spools  designated  as 
008  and  019.  Each  composite  strand  consisted  of  3000  fibers  with  a 
nominal  cross-sectional  area  of  .66  mm^.  Samples  with  different  gauge 
lengths  were  tested,  the  majority  of  which  were  2  and  10  inches  for  spool 
008  samples,  and  2  inches  for  019  spool  samples.  Ten  inch  samples  from 
the  019  spool  were  tested  in  previous  investigations  and  are  also  included 
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for  interpretation..  Two  additional  samples,  one  5  inch  and  one  16  inch 
gauge  length  were  also  tested  and  this  data  used  for  compliance  calculations 
in  addition  to  normal  data  reduction. 
2.     Test  Method  Development 

The  samples  were  required  to  be  made  such  that  the  strand  could  be 
gripped  and  held  under  load  in  the  testing  instrument.  In  addition,  it  was 
desired  that  the  strand  fail  in  the  gauge  length  area  and  not  in  the  gripping 
area.  Two  different  testing  methods  were  identified  and  investigated  to 
accomplish  these  requirements.  The  first  method  involved  securing  the  test 
strand  ends  in  a  copper  tube  (.0669  inch  OD  and  .0390  inch  ID)  utilizing 
an  epoxy  type  adhesive  with  the  tubes  separated  by  the  desired  gauge 
length.  The  two  ends  were  then  be  placed  in  an  adapter  collet  (from  a 
manufactured  jewler's  lath  having  .0700  inch  ID  and  .500  inch  OD)  which 
in  turn  were  placed  in  the  grips  of  the  testing  instrument.  Several 
problems  were  encounter  in  the  implementation  of  this  method.  First,  the 
adapter  collet/copper  tube  interface  area  was  too  small,  producing  a  high 
stress  concentration  area  on  the  strand  itself  upon  gripping,  causing  failure 
of  the  strand  in  the  grip.  Also,  the  adapter  collets  were  not  designed  to 
accept  the  type  of  loading  the  gripping  procedure  produced  and  they  failed 
after  a  relatively  low  number  of  loadings.  The  second  method  involved 
using  one-half  inch  aluminum  bar  stock  to  make  1/2  inch  long  pellets  that 
were  secured  to  the  ends  of  the  strands  with  the  epoxy,  resulting  in  a 
"dumbbell"  type  configuration  of  the  test  strand.  Each  pellet  had  a  #60 
drill  size  hole  (dimensions)  machined  into  the  center  of  it  to  accommodate 
the  ends  of  the  strands.  This  configuration  allowed  the  test  samples  to  be 
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placed  directly  into  the  testing  instrument,  thereby  simplifying  the  test 
procedures.   Jigs,  which  were  made  out  of  aluminum  angle  iron,  were 
designed  and  manufactured  to  maintain  the  pellets  at  the  proper  gauge 
length  during  the  adhesive  bounding  process.  At  the  same  time,  the  jigs 
protected  the  samples  from  handling  damage.  This  was  the  procedure  that 
was  followed  to  make  the  test  articles  and  is  explained  in  further  detail  in 
the  following  sections. 

Application  of  hydrostatic  pressure,  via  the  hydraulic  grips,  to  the 
pellets  produced  a  hydrostatic  state  of  compressive  stress  in  the 
sample/aluminim  pellet  interface.  This  increased  the  shear  strength  of  that 
portion  of  the  strand  which  faciliated  the  transfer  of  the  external  load  to 
the  fiber  filaments  within  the  free  strand.  Through  trial  and  error,  it  was 
determined  that  a  reduction  of  .0025  to  .0035  inches  of  the  outer  diameter 
of  the  aluminum  pellet  upon  application  of  gripping  pressure  was  sufficient 
to  produce  an  effective  hydrostatic  pressure.   This  amount  of  deformation 
was  also  the  limit  to  prevent  permanent  deformation  of  the  collets  inside 
the  grips  themselves.  To  consistently  obtain  the  desired  yield 
characteristics,  the  pellets  were  all  annealed  at  450  deg  Celsius  for  1  1/2 
hours  and  then  oven  cooled  overnight.  This  not  only  produced  the  desired 
deformation  properties  but  also  acted  to  rid  the  pellets  of  any  residual  oil 
and  dirt  deposited  during  the  machining  process.  The  annealed  pellets 
were  then  tested  by  applying  incremental  steps  of  hydraulic  pressure  to  the 
grips  holding  the  pellet  to  determine  the  gauge  pressure  in  the  testing 
system  that  produced  the  desired  deformation.  This  optimal  pressure  was 
estimated  to  be  500  psi  indicated. 
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3.   Stress  Wave  Attenuation 

The  intended  purpose  of  a  stress  wave  attenuator  was  to  prevent 
secondary  failures  of  the  strand  sample  caused  by  the  stress  wave 
accompanied  by  the  initial  failure.    Immediately  after  the  initial  failure  of 
the  sample  the  tensile  energy  stored  in  the  unbroken  portion  of  the  sample 
is  released  resulting  in  a  tensile  stress  wave  traveling  toward  the  two  ends 
of  the  grip.  Upon  reflection  from  the  strand/pellet  grip  interface,  the 
tensile  wave  is  converted  to  a  compressive  wave.  Since  the  strand  sample  is 
long  and  slender,  the  compressive  wave  almost  always  caused  a  secondary 
failure  at  the  strand/pellet  interface  due  to  compress  buckling.  Additional 
secondary  failure  sites  throughout  the  length  of  the  sample  were  also 
fequently  observed,  therby  precluded  a  definitive  identification  of  the 
original  initial  failure  site.  It  was  experimentally  observed  that  when  the 
tensile  stress  was  attenuated  during  the  intial  propagation  stage,  the 
secondary  breakages  could  be  eliminated.  Because  in  the  current  testing 
configuration,  the  test  sample  was  gripped  in  a  vertical  position,  sand  and 
oil  were  ruled  out  as  candidates  as  stress  wave  attenuators,  even  though 
both  of  these  have  proven  successful  in  the  past.1  The  use  of  paraffin  wax 
was  also  attempted  but  without  success,  as  was  modeling  clay.  The  material 
adopted  for  these  tests  was  a  commercially  available  generic  hot-melt 
adhesive  that  was  applied,  via  a  hot  melt  glue  gun,  to  the  entire  gauge 
length  of  the  test  specimen.  The  general  chemical  constituents  of  the 


1    Wu/Nypiuk    testing    at    Livermore   Lab 
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adhesive  include  a  styrene  derivative,  a  hydrocarbon  resin,  artificial  wax, 
and  ethylene  vinyl  acetate  coploymer  resins. 
4.   Test  Specimen  Preparation 

The  composite  strands  were  originally  manufactured  in 
approximately  18  inch  sections  and  cut  to  lengths  that  were  1  1/4  inches 
longer  that  the  desired  gauge  length  by  the  use  of  a  Dremel  tool  (using  a 
cutoff  grinding  wheel  attachment).  The  additional  length  was  used  to 
insure  the  strand  extended  through  the  entire  length  of  the  pellet.  Once  the 
required  number  of  samples  were  cut  (8  samples  per  batch),  each  end  was 
cleaned  by  rinsing  them  through  two  separate  acetone  baths  using  a  glass 
syringe.   Additionally,  the  utensils  (stirring  stick  and  the  syringe  used  for 
the  placement  of  the  epoxy)  were  also  cleaned  in  the  same  manner.  An 
epoxy  that  consisted  of  55%  by  weight  Dow  DER-332  Epoxy  Resin  and 
45%  by  weight  Texaco  Jeffamine  T403  was  then  made  and  stirred  for  a 
minimum  of  3  minutes  to  ensure  thorough  mixing.  The  epoxy  was  then 
placed  into  a  vacuum  chamber  to  rid  it  of  the  air  bubbles  developed  while 
mixing.  This  epoxy  was  used  because  of  its  particularly  good  wetting 
properties  and  its  slow  setting  time.  This  allowed  sample  preparation  to  be 
conducted  in  a  slow,  methodical  process  reducing  the  possibility  of 
damaging  the  samples  during  manufacturing. 

The  pellets  that  were  attached  to  the  end  of  the  strands  (2  each) 
were  then  prepared  by  placing  cellophane  tape  on  three  of  the  four  ends  of 
each  set  of  pellets  with  small  holes  punctured  through  the  tape  in  the  center 
of  each  machined  hole  in  the  pellet.  The  use  of  a  pin  and  microscope 
helped  in  placing  the  puncture  hole  in  the  desired  location.  This  process 
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was  done  to  help  center  the  strand  in  the  pellet  during  the  production  of  the 
test  samples.  The  pellets  were  then  placed  in  the  jigs  with  the  pellet  with 
tape  on  both  sides  being  the  bottom  pellet  and  the  top  pellet  having  the 
taped  side  down.  They  were  held  in  proper  place  in  the  jig  by  rubber 
bands,  allowing  for  the  easy  removal  (cutting  of  the  rubber  bands)  of  the 
test  sample  from  the  jig  at  time  of  testing.  The  rear  portion  of  the  angle 
iron  type  jigs  was  machined  away  enabling  access  to  the  rear  of  the  strand, 
making  the  application  of  the  hot-melt  adhesive  easier.  Once  the  jigs  and 
the  strands  were  prepared  and  ready,  the  epoxy  was  placed  into  a  plastic, 
reusable  syringe  with  a  modified  tip  consisting  of  the  bottom  tip  of  a  5  mm 
mechanical  pencil.  This  modification  allowed  the  injection  of  epoxy 
directly  into  the  top  holes  of  the  pellets  until  the  epoxy  flowed  out  the 
bottom.  Since  test  samples  were  made  in  batches  of  eight,  each  of  the  16 
pellets  were  filled  with  epoxy  in  this  manner.  The  strands  were  then 
inserted  into  the  pellets  such  that  the  strand  extended  through  the  entire 
length  of  each  pellet.  The  jig  was  then  placed  in  a  vertical  stand.  Surface 
tension  retaining  the  liquid  epoxy  in  the  cavity  of  the  pellet  alleviated  the 
problem  of  inadvertently  attaching  the  pellet  to  the  jig.  Once  all  eight 
samples  were  made,  the  stand  was  placed  in  an  oven  which  was  pre-heated 
to  58  deg  Celsius  (to  accelerate  cross-linking  of  the  epoxy)  and  cured  for 
24  hours.  After  curing  of  the  adhesive  in  the  pellet,  each  sample  had  the 
hot-melt  adhesive  applied,  by  hand,  to  the  entire  length  of  the  strand.  The 
samples  were  then  ready  for  testing.  Figure  8,  Appendix  A  shows  the 
completed  samples. 
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C.    TESTING  PROCEDURES 

1.  Alignment 

Before  any  testing  was  conducted,  the  test  instrument's  upper  and 
lower  grips  were  aligned  by  griping  a  precision  ground  1/2  inch  diameter 
steel  drill  rod  in  the  grips  with  the  grip  surfaces  separated  by 
approximately  1/4  inch.  A  load  of  1000  Kg  (the  load  cell  design  limit)  was 
then  carefully  applied  and  maintained  (ensuring  the  load  never  exceeded 
design  limits).  To  eliminate  all  of  the  mechanical  play  or  looseness  in  the 
system,  all  of  the  rigid  coupling  attachments  were  then  tightened  under  this 
tension,thereby  setting  the  alignment  of  the  grips.  To  ensure  a  consistent 
mechanical  reference  state  for  every  sample  tested,  the  grip  alignment  was 
checked  and  verified  approximately  one-half  way  through  each  testing 
interval. 

2.  Calibration 

Load  cell  calibration  was  conducted  before,  during  and  after  the 
testing  period.   Calibration  was  conducted  using  Instron  calibration  weights 
(traceable  to  National  Bureau  of  Standards)  with  load  increments  of  5 
Kg.   Calibration  data  is  presented  inTable  I,  Appendix  B.  In  addition,  a 
time  drift  or  stability  check  of  the  system  was  conducted  over  a  24  hour 
period  utilizing  a  constant  load  of  15  Kg.  The  greatest  drift  recorded 
during  the  period  was  within  .17  per  cent  of  the  original  load.   As  shown, 
the  calibration  remained  well  within  the  published  specification  limit  of 
±1.0%. 
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3.   Testing 

Testing  of  the  prepared  samples  involved  transferring  the  samples 
from  the  jigs  to  the  test  instrument  and  initializing  and  starting  the  test  via 
the  computer.  The  fully  automated  testing  system  recorded  all  the  data  and 
stopped  the  test  upon  failure  of  the  sample. 

The  placement  of  the  sample  into  the  test  instrument  required  two 
people,  one  to  assist  the  pellet  of  the  strand  in  sliding  into  the  collets 
internal  to  the  grips  and  the  other  to  operate  the  console  that  controls  the 
crosshead  displacement.  The  sample  was  first  placed  into  the  upper  grip 
and  inserted  far  enough  into  the  collet  such  that  the  lip  of  the  pellet  was 
inserted  beyond  the  lip  of  the  collet.  This  ensured  that  the  gripping  force 
was  distributed  over  the  entire  surface  area  of  the  pellet.   A  gripping  force, 
measured  by  a  hydraulic  gauge  pressure,  of  approximately  500  psi  was 
then  applied  to  the  grip  housing  via  the  hydraulic  pump.  The  crosshead 
was  next  lowered  (manual  mode,  fast  speed)  such  that  the  bottom  pellet  was 
approximately  1/4  inch  above  the  upper  surface  of  the  lower  grip  housing. 
At  this  point,  the  speed  select  was  changed  to  manual  mode,  slow  speed  and 
a  load  balance  was  conducted.  Then  the  bottom  pellet  was  inserted  into  the 
collet  of  the  lower  grip  to  the  same  relative  position  as  previously 
described  and  gripping  pressure  applied.    The  crosshead  was  adjusted  such 
that  a  slight  negative  load  (0.0  to  -.10  Kg)  was  indicated.  The  crosshead 
position  reference  was  then  reset  to  zero  completing  the  installation. 

The  computer  was  initialized,  dependent  on  gauge  length,  utilizing 
the  appropriate  test  method  program.   Sample  identification  numbers  were 
developed  such  that  they  could  be  easily  identified.  For  example,  sample 
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01902-15  was  the  15m  sample  of  2  inch  gauge  length  of  the  AS4-019  class 
of  composite  strand  that  was  tested.  Once  the  required  information  of  the 
interactive  program  was  entered,  the  test  was  conducted.  After  testing  to 
failure,  the  sample  was  removed  from  the  test  instrument  and  examined  to 
determine  the  location  of  failure.  A  summary  of  all  pertinent  data 
collected  during  each  test  was  then  documented  in  the  computer  generated 
test  report.    A  sample  test  report  is  presented  in  Figure  9,  Appendix  A. 
The  raw  data  was  finally  stored  on  the  hard  disk  and  backed  up  on  a  floppy 
disk. 
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IV.      RESULTS 

A.  TEST  DATA 

1.  Current  Testing 

A  total  of  82  strands  were  tested  to  failure  of  which  56  were  from  the 
AS4-008  graphite  spool  and  26  from  the  AS4-019  spool.  Load  and 
displacement  at  failure  data  for  all  tests  are  presented  in  Tables  II  thru  VI, 
Appendix  B.  Compliance  calculation  methodology  for  the  correction 
factor  applied  to  sample  displacements  due  to  system  mechanical 
displacement  is  discussed  in  Appendix  C. 

2.  Previous   Testing 

Data  from  previous  tests  of  008  and  019  composite  strands  that 
were  part  of  the  same  physical  set  of  the  current  test  samples  has  been 
included  in  the  data  base  for  use  in  the  analysis  of  data  in  this  report.  This 
was  done  to  better  define  the  reliability  curve  that  the  process  of  data 
reduction  produces.  The  source  of  this  data  is  from  tests  conducted  at 
Lawrence  Livermore  National  Laboratories  by  Dr.  Edward  M.  Wu  and 
Mr.  Glenn  Nypiuk  [Ref.  6].  This  data  is  presented  in  Tables  VII  and  VIII, 
Appendix  B.  AS4-008  data  from  these  tests  are  designated  set  A  as 
compared  to  AS4-008  data  set  B  which  designates  current  test  data. 

Fiber  strength  research  for  both  (008  and  019)  graphite  spools  was 
conducted  by  Dr.  Edward  M.  Wu  and  Mr.  Nypiuk  [Ref.  6],  Lt.  David 
Keith  Bell,  USN  [Ref.  7],  LCDR  Carl  Engelbert,  USN  [Ref.  8]  and  Mr.  Jim 
Nageotte  of  the  Advanced  Composites  Laboratory,  Naval  Postgraduate 
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School.  The  fiber  data  from  the  above  research,  which  is  to  be  used  in  the 
final  analysis  of  the  present  research,  is  contained  in  Ref.  9. 

B.  HOT-MELT  ADHESIVE  CONTRIBUTION  TO  SAMPLE 
STRENGTH 

By  using  hot-melt  adhesive  as  a  stress  wave  attenuator,  the  question 

arises  that,  what  contribution,  if  any,  does  the  hot-melt  adhesive  make  to 

the  strength  of  the  sample?  This  issue  was  resolved  by  testing  the  adhesive 

and  determining  its  modulus  which  was  then  compared  to  that  of  the 

strand.  Graphical  representation  of  the  hot-melt  load-displacement  data  is 

presented  in  Figure  10,  Appendix  A.   Knowing  the  cross-sectional  area  of 

the  test  sample,  the  modulus  of  the  adhesive  was  calculated  to  be  2.845E3 

lbs  /  in^.  When  compared  to  the  modulus  (approximated  at  20E6  lbs  /  in^) 

of  the  composite  strand,  this  is  less  than  .014  per  cent  of  the  strand 

modulus.  In  addition  to  the  modulus  comparison,  experiments  were 

conducted  to  determine  the  effect  the  addition  of  the  adhesive  to  the  strand 

would  make  on  the  raw  data.  This  was  done  by  first  testing  a  single  strand 

(without  adhesive)  to  a  load  of  approximately  50%  of  the  failure  load. 

This  load  application  procedure  was  repeated  several  times  to  ensure  that 

all  the  fibers  that  were  weak  and  would  break  under  a  load  up  to  and 

including  the  peak  applied  load,  actually  broke.  This  resulted  in  a  fixed, 

repeatable,  load-displacement  curve  for  the  test  sample,  indiscriminate  of 

how  many  times  a  load  (of  the  same  peak  magnitude)  was  applied.  The 

hot-melt  adhesive  was  then  applied  to  the  strand  and  the  sample  tested  to 

the  same  load  as  before.  As  can  be  seen  in  Figure  1 1 ,  Appendix  A,  the 
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load-deformation  curve  of  the  sample  with  the  adhesive  lies  directly  on  top 
of  the  same  curve  of  the  sample  without  the  adhesive. 

In  analyzing  the  results  of  both  of  these  experiments,  it  was 
concluded  that  in  the  context  of  the  testing  done  to  the  composite  strands  on 
this  project,  that  the  hot-melt  adhesive  did  not  add  to  the  strength  of  the 
composite  strand  and  therefore  the  use  of  the  hot-melt  adhesive  as  a  stress 
wave  deflector  did  not  effect  the  data. 

C.    DATA   ANALYSIS 

1.   Data  With  Experimental  Artifact 

In  this  type  of  data  analysis,  the  data  of  those  strands  that  failed  in 
the  gauge  length  portion  of  a  test  sample  are  considered  the  intrinsic  load 
of  that  particular  sample.  The  question  then  arises  of  how  should  the  data 
of  the  samples  that  have  multiple  failure  sites  including  one  failure  site  at 
the  grip/strand  interface,  be  handled?  We  classified  this  data  as  censored 
data  (due  to  the  potential  experimental  artifacts  which  may  have  lowered 
the  strength)  .  Based  on  this  physical  consideration,  the  failure  load  (lJ 
for  censored  data  samples  must  be  always  less  than  and  at  best  equal  to  its 
intrinsic  value,  that  is,  (Lj)  or  (lJ  <  (Lj)  .  Our  goal  is  to  be  able  to  use  this 
censored  data  in  combination  with  the  intrinsic  data  in  the  data  analysis. 
To  achieve  this,  during  data  analysis  a  procedure  using  the  Maximum 
Likelihood  Estimator  (MLE)  was  utilized.  The  MLE  concept  for  a  two 
parameter  Weibull  distribution  is  discussed  in  Appendix  E. 

In  data  analysis,  it  is  desired  to  plot  the  Weibull  probability  of 
failure  of  the  material  F*(Lj)  versus  its  intrinsic  load.  To  do  this,  one 

must  order  the  experimental  data  (failure  load)  in  ascending  order  and 
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assign  a  rank  to  each  point.  In  this  analysis,  expected  ranking  was  used 
such  that  the  rank  of  (Lj)  is 

forx=  1,2,3 N 


N+l 

where  x  is  the  numerical  sequence  of  the  failure  load  after  ordering  and  N 
is  the  total  number  of  samples  tested.  Because  some  of  the  data  is  censored 
and  not  intrinsic,  a  MLE  analysis  is  conducted  to  estimate  the  model  of,  or 
the  parameters  of,  the  Weibull  distribution  associated  with  this  particular 
set  of  data.  A  new  probability  of  failure  F  [L-)  is  then  determined  for  each 
point  of  censored  data  (lJ  as  shown  in  Figure  12,  Appendix  A.   In 
analyzing  F (L-),  we  note  that  if  F  (lJ >  F (lj,  it  physically  implies  that  the 
failure  load  (lj  was  greater  than  its  intrinsic  load  (L;)  or  (lJ  >  (L;),  which 
is  inadmissible  in  the  physical  model.  Therefore,  in  this  case,  for  each 
point  where  F  (lJ  >  F  [L-J,  we  assume  (lj  =  (L;)  and  F  (lJ  =  F(Lj)  or  that 
the  failure  load  is  the  intrinsic  load.   If  F  (lJ  <  F  (lJ,  the  value  of  F  (lJ  is 
carried  forward  in  the  calculations.   The  new  values  of  F^Cj  are  then 
combined  with  the  original  ranking  of  the  intrinsic  loads  and  the  rank  (with 
rank  remaining  keyed  to  the  load)  is  reordered  in  ascending  order.  It 
should  be  noted  that  the  MLE  process  should  be  performed  on  a  sample  set 
of  a  single  dimension  only.  This  is  based  on  the  fact  that  when  the  range  of 
the  probability  of  failure,  F(Lj)  {ie,  the  sample  size  N}  is  small,  the 
assumption  of  a  linear  Weibull  distribution  is  a  good  approximation  of  the 
model.  In  other  words,  the  model  is  assumed  to  be  piecewise  linear  when 
the  sample  size  is  small. 
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2.   Data  Considered  Intrinsic 

This  type  of  analysis  assumes  that  the  failure  load  of  a  particular 
test  sample  is  the  intrinsic  load  of  the  sample,  or  the  sample  failed  at  its 
statistically  distributed  intrinsic  stress  level  or  load.  This  assumption 
excludes  any  experimental  artifact  that  may  be  associated  with  the  data. 
This  is  the  basis  of  the  analysis  used  in  this  investigation. 

a.    Data   Interpretation 

Since  both  composite  and  fiber  strength  are  functions  of  gauge 
length,  a  mathematical  model  must  be  determined  to  normalize  all  the  data 
to  one  single  metric  size  so  that  a  direct  comparison  of  all  the  data  can  be 
made.  The  standardization  by  a  metric  dimension  can  be  chosen  arbitrarily 
without  affecting  the  relative  relations.  The  metric  selected  by  the  current 
data  analysis  is  one  centimeter.  The  method  for  this  standardization  to  a 
common  metric  dimension  is  discussed  in  Appendix  D.  The 
standardization  used  in  the  data  analysis  of  this  study  is  the  vertical  shift 
type. 

To  determine  if  the  data  of  like  materials  and  gauge  lengths 
are  of  the  same  set  (identical  test  samples),  the  Weibull  probability  of 
failure  plots  are  compared  for  each  batch  of  tests  completed.  Figures  13, 
14,  and  15,  Appendix  A  show  a  comparison  of  each  batch  of  data  collected. 
In  the  analysis,  batch  one  was  the  first  half  of  the  sample  set  which  was 
compared  to  the  second  half  of  the  sample  set.  Batch  two  compared  all  the 
odd  number  test  samples  to  that  of  the  even  number  test  samples.  The  data 
band  width  is  dependent  on  the  experimental  circumstances.  A  narrow 
band  implies  that  the  experimental  techniques  among  the  batches  are 
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similar,  justifying  the  merging  of  the  different  test  batches  into  one  data 
set.  As  shown,  the  data  band  width  for  each  of  the  sample  sets  is  narrow, 
indicating  that  the  data  belongs  to  the  same  set,  which  is  consistent  to  the 
actual  experimental  circumstances.  By  confirming  that  like  batches  are  the 
same,  they  can  be  merged  to  produce  single  sample  sets  of  like  materials 
and  gauge  lengths.  Figures  16  through  20,  Appendix  A  presents  the 
merged  data  of  the  sample  sets  as  compared  with  the  assumed  Weibull 
distribution  model  obtained  with  the  use  of  a  Maximum  Likelihood 
Estimator  as  described  in  Appendix  E.  The  band  width  of  the  two  sets  of 
008  10"  data  (sets  A  and  B)  was  then  analyzed  as  shown  in  Figure  21, 
Appendix  A  and  the  two  sets  determined  to  be  of  the  same  set.  The  008 
10"  data  was  then  merged  and  compared  to  the  assumed  Weibull 
distribution  model  as  shown  in  Figure  22,  Appendix  A.  The  values  of  a 
and  (3  produced  by  the  MLE  method  for  each  set  of  data  are  presented  in 
Table  IX,  Appendix  B.  All  the  data  was  then  converted  to  a  1  cm  metric 
format  to  enable  a  direct  comparison  with  fiber  data  of  the  same  metric. 
Figures  23  and  24,  Appendix  A  present  this  data  in  terms  of  composite 
sample  strength  (per  bundle)  and  Figures  25  and  26,  Appendix  A  present  it 
in  terms  of  composite  strength  per  fiber  in  a  bundle.   As  shown  in  Figure 
27,  Appendix  A,  in  a  direct  comparison  of  composite  strength  (per  fiber), 
the  composite  made  with  the  008  spool  of  AS4  fiber  appears  to  be  stronger 
than  that  made  with  the  019  spool  of  fiber  even  though  both  spools  came 
from  the  same  manufacturing  line  and  were  assumed  to  be  identically 
produced.  This  is  an  indication  that  fiber  spools  that  may  have  even  been 
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produced  in  the  same  line,  let  alone  the  same  lot,  have  to  be  regarded  as 
having  different  statistical  characteristics  and  treated  as  such. 

Much  the  same  type  of  analysis  can  be  done  with  the  fiber 
data.  Figures  28  and  29,  Appendix  A  show  the  band  width  of  the  fiber 
data,  which  confirms  that  the  data  belong  to  the  same  set  of  like  materials 
and  gauge  lengths.  Figures  30  and  31,  Appendix  A  present  the  merged 
fiber  data  as  compared  with  the  assumed  Weibull  distribution  model.  The 
values  of  the  trimodal  Weibull  distribution  for  this  data  are  presented  in 
Table  IX,  Appendix  B.   As  noted  in  Figure  30,  Appendix  A,  the  019  fiber 
probability  of  failure  is  judged  to  be  bimodal,  or  the  lower  tail  having 
different  Weibull  distribution  parameters  than  the  upper  tail.  The  modal 
shape  of  the  data  was  determined  by  visual  examination  of  merged 
individual  ordered  data  plotted  on  a  Weibull  axis.  The  data  was  then 
identified,  in  this  case,  as  either  unimodal  or  bimodal.  The  fitting  of  the 
data  to  the  upper  and  lower  curves,  including  the  identification  of  the 
lower  transition  point  was  performed  heuristically.   The  MLE  program 
was  used  to  determine  the  values  of  a  and  p  for  both  the  lower  and  upper 
curves.  The  use  of  MLE  is  permitted  for  the  bimodal  case  because  it  is 
considered  to  be  a  piecewise  linear  distribution  in  the  Weibull  plot  with  a 
distinct  transition  point.       Figures  32  and  33,  Appendix  A  present  the 
Weibull  probability  of  failure  of  the  fiber  when  converted  to  a  metric  of  1 
cm.  Figure  34,  Appendix  A  shows  the  comparison  of  the  fiber  strengths 
and  is  consistent  with  Figure  27,  Appendix  A. 
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b.    Local  Load  Sharing  Analysis 

The  data  from  Figures  25  and  32,  Appendix  A  were  combined 
to  produce  Figure  35,  Appendix  A.  Likewise,  Figures  26  and  33  were 
combined  to  produce  Figure  36,  Appendix  A.  These  plots  are  used  in 
conjunction  with  the  Local  Load  Sharing  (LLS)  Prediction  Model  with 
extensions,  to  enable  prediction  of  composite  reliability  from  trimodal 
Weibull  fiber  statistics  [Ref.  10].  The  model  requires  two  inputs,  the  fiber 
model  parameters  and  the  ineffective  length,  in  order  to  predict  the  model 
for  the  composite  reliability.   The  ineffective  length,  however,  is  not 
known,  so  the  LLS  program  [Ref.  10]   was  written  in  a  normalized  (Xi/p) 
format.   Therefore,  when  the  model  parameters  of  a  particular  set  of  fiber 
data  are  used  in  the  program,  the  LLS  curve  that  is  produced  is  referenced 
to  the  normalized  composite  data  of  the  same  gauge  length.  In  this 
analysis,  one  set  of  data  (ie,  the  008  data)  was  used  in  the  model,  since  both 
the  fiber  and  composite  statistics  (Figure  36,  Appendix  A)  were  known. 
Fiber  statistics  were  entered  using  the  model  parameters  in  a  trimodal 
Weibull  distribution  format.  The  model  then  generates  a  predicted  value 
for  the  composite  reliability  from  which  the  ineffective  length  is  estimated. 
This  is  done  by  shifting  the  predicted  reliability  model  horizontally  to  the 
right  until  the  model  approaches  the  value  of  the  associated  composite  test 
data.  Since  the  amount  shifted  is  a  function  of  the  ineffective  length  of  the 
composite,  it  is  assumed  to  be  the  same  for  a  composite  made  of  the  same 
fiber  (AS4)  but  of  different  spools.  The  same  magnitude  of  shift  is  then 
used  in  conjunction  with  the  019  fiber  statistics  and  the  LLS  model  to 
predict  the  019  composite  reliability.  This  prediction  can  be  compared  and 
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cross-plotted  directly  with  Figure  35,  Appendix  A,  making  analysis  of  the 
LLS  Prediction  Model  possible.  The  process  can  also  be  reversed  by  using 
the  008  data  to  obtain  for  the  ineffective  length  and  then  predicting  the  019 
composite  reliability  with  the  LLS  model,  which  can  then  be  compared 
with  Figure  35,  Appendix  A.  Figures  37  and  38,  Appendix  A  presents  the 
graphical  comparison  of  the  predicted  composite  reliability  and  the 
experimental  test  data  for  both  008  and  019  sample  sets. 

As  can  be  seen  in  Figures  37  and  38,  the  LLS  composite 
reliability  prediction  exhibits  outstanding  correlation  with  the  experimental 
data.  In  addition,  it  is  sufficiently  sensitive  to  predict  the  shifted  composite 
strength  due  to  the  abnormally  weak  fiber  tail  associated  with  the  019 
fiber.  It  can  also  be  inferred  that  the  LLS  reliability  prediction  provides 
the  upper  bound  on  reliability  since  the  model  assumes  a  well-made 
composite  with  a  consistent  ineffective  length.  This  is  further  supported  by 
the  fact  that  the  experimental  data  of  the  test  specimens  (which  are 
considered  well-made)  correlated  so  closely  with  the  predicted  model. 

D.    SUMMARY 

In  summary,  the  major  effort  of  this  investigation  was  in  the  collection 
of  a  comprehensive  set  of  composite  data  to  be  interpreted  in  conjunction 
with  parent  fiber  data  sets  previously  collected.  The  model  for  the  data 
interpretation  is  the  Local  Load  Sharing  model  which  uses  fiber  statistics  to 
predict  the  spatial  density  and  clustering  of  the  fiber  failure  sites.  From 
the  local  stress  concentrations  that  arise  from  load  sharing,  the  probability 
of  composite  failure  is  predicted.  This  analytical  prediction  is  then 
compared  to  the  composite  strength  data  sets  obtained  in  this  investigation. 
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Since  the  constituent  fiber  statistics  are  the  benchmark  for  composite 
reliability  prediction  using  the  LLS  prediction  model,  it  is  important  to 
ensure  that  the  experimental  fiber  data  was  obtained  in  such  a  manner  that 
the  data  accurately  reflects  the  statistical  distribution  of  the  fiber  strength. 
This  was  accomplished  in  previous  investigations  by  testing  a  sufficiently 
large  sample  population  which  ensures  that  the  modal  occurrence  and  the 
trimodal  Weibull  distribution  parameters  estimated  are  representative  of 
the  underlying  fiber  population.  In  addition,  the  use  of  multiple  operators 
in  conducting  the  tests  also  prevented  the  inadvertent  skewing  of  the  data. 
Finally  different  spools  of  fiber  were  used  to  identify  normal 
manufacturing  variables  that  are  always  present  during  the  production 
process. 

To  verify  the  prediction  model  with  composite  test  data,  the  composite 
data  must  again  be  beyond  question  in  proper  representation  of  composite 
strength  statistics.  In  this  investigation,  much  time  was  spent  in  developing 
the  test  method  procedures  to  ensure  that  the  experimental  data  was 
accurate.  The  test  sample  gripping  technique  developed  and  utilized  for  the 
longitudinal  tension  testing  was  designed  to  prevent  the  normally  dominant 
shear  failure  mode  associated  with  tension  testing  from  causing  the  test 
specimen  failure.  This  was  accomplished  with  a  unique  encapsulating  grip 
that  applyed  hydrostatic  pressure  to  the  sample  at  the  grip/sample  interface 
which  was  induced  by  static  compression.  The  use  of  the  hot-melt  adhesive 
to  attenuate  the  stress  wave  that  accompanies  the  sample  failure,  prevented 
the  shattering  of  the  sample  and  enabled  the  identification  of  the  failure 
characteristics.  In  addition,  different  sample  dimensions  (gauge  lengths) 
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were  tested  to  extend  the  data  range  and  enable  the  use  of  the  methodology 
of  data  interpretation  by  shifting  the  weakest  link.  All  this  was  successfully 
accomplished,  producing  great  assurance  in  the  experimental  results. 

In  interpretation  of  the  data,  for  the  case  where  the  fiber  is 
unimodal,  the  parameters  for  the  two  parameter  Weibull  distribution  were 
estimated  using  the  maximum  likelihood  estimator.  For  the  case  of  data 
(019  spool  of  fiber)  which  displayed  a  bimodal  distribution  a  model 
consisting  of  three  piecewise  linear  Weibull  distributions  was  used.  The 
parameters  for  the  trimodal  Weibull  distribution  model  were  estimated 
from  the  ordered  data  set  with  censoring  to  the  left  and  right  respectively. 
The  composite  test  samples  that  were  used  in  this  investigation  were 
considered  to  be  exceedingly  well-made  composite  strands  with  uniform 
fiber  and  matrix  adhesion  and  therefore,  a  uniform  ineffective  length.  The 
ineffective  length  is  controlled  by  the  interfacial  strength  (a  function  of  the 
manufacturing  process)  and  the  shear  modulus  of  the  matrix  (materials). 
Composites  with  irregular  fiber-matrix  adhesion  and  voids  have  large 
ineffective  lengths,  the  limiting  case  being  a  composite  with  no  matrix  (a 
bundle).  It  was  observed  that  for  this  set  of  well-made  composite,  the 
Local  Load  Sharing  Model  of  reliability  prediction  exhibited  outstanding 
correlation  with  experimental  data  and  was  sufficiently  sensitive  to  predict 
deficient  composite  strength  due  to  a  specific  fiber  spool  with  an 
abnormally  weak  lower  tail.   In  addition,  it  provided  an  upper  bound  on 
the  composite  reliability.  The  documented  weaker  lower  tail  of  the  019 
fiber  was  predicted  in  the  composite  reliability  by  the  LLS  model, 
verifying  six  orders  of  extrapolation  of  size  in  reliability  prediction.  It 
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should  be  noted  that  all  analytical  results  were  consistent  with  experimental 
observation,  not  only  in  trend  but  also  in  substantial  agreement  with 


magnitudes. 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  investigation  suggest  that  the  Local  Load  Sharing 
composite  reliability  prediction  method  using  a  trimodal  Weibull 
distribution  is  a  viable  method  of  predicting  composite  structural  reliability 
using  fiber  strength  statistics.  In  design,  this  methodology  provides  a 
relation  of  strength  data,  based  on  small  laboratory  sample  dimensions,  to 
actual  large  structural  dimensions  for  a  fail  safe  design.  In  materials 
development,  this  methodology  can  be  used  to  identify  the  parametric  role 
of  fiber  and  matrix  strength  properties  for  possible  improvements.   In 
manufacturing,  this  methodology  can  be  used  to  prescreen  fiber  for  zero- 
reject  fabrications.   In  maintenance  and  repair,  this  methodology  provides 
quantitative  guidelines. 

The  following  areas  were  not  fully  treated  in  this  investigation  and  are 
recommendations  of  follow-on  research: 

1 .  That  larger  physical  dimension  test  samples  be  tested  such  that 
next  several  orders  of  data  extrapolation  can  be  verified. 

2.  That  a  mathematical  formalism  of  the  statistical  methods  to 
estimate  the  parameters  for  the  trimodal  Weibull  distribution  be 
established. 

3.  That  a  mathematical  formalism  of  the  statistical  methods  to 
estimate  the  best  fit  ineffective  length  based  on  the  composite  and  fiber  data 
be  conducted. 
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4.     That  a  closed  form  approximation  of  the  Local  Load  Sharing 
composite  reliability  prediction  method  using  trimodal  Weibull  distribution 
parameters  be  developed. 
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APPENDIX  A 
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Figure  1 .  Composite  Failure  Model 
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X  -  Load  Per  Fiber 

d  -  Fiber  Element  Length 

Kr  -  Load  Concentration  Factor 

N  -  Number  Of  Fibers 

M  -  Number  Of  Bundles 


-  Fiber  Break 


Figure  2.  Local  Load  Sharing  Model 
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Figure  3.   Instron  Testing  Instrument 
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Figure  5.  Test  Method  -  5"  Gauge  Length 
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Figure  6.  Test  Method  -  10"  Gauge  Length 
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!  Reverse  load  polarity,"  no 

!  Sailing  rate:    6.67  pt./sec, 

!  disable  reduct 

!  Load-cell  used:   1000.00  Kgrans 
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Data  GOHpression:  disabled  ! 
Machine  control  parameters:  

disable  lest  stop  at: 
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Figure  7.  Test  Method  -  16"  Gauge  Length 
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Figure  8.  Composite  Strand  Test  Samples 
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Naval  Fostoradi'ate  Schcpl 
Composites  Mater i a  1  Lab 
Thesis  -  Johnson 

Sample  Identification:  00916-01  Test  Date:  Apr i 1  II, 1991 

Cut  of  I  specimens,  0  excluded. 

Sample  consents:  Very  good  break,  ends  intact,  rest  shattered. 

Displcifent  Load     Stress     ?'r?jn  X  Strain  Displc«?nt  Load     Stress  Strain 

at  at        at        at  at  at  Auto.  at  Auto.    at  Auto.  at  Auto. 

?p?ci«ien    Mar. Load  Mar. Load  Mar. Load  Ha^:  .Lead  Mar. Load  Break  Break     Break     Freak 

Number       (»i)  (Kg)  (Kg/M««2>  lim/*-.)  (X)  tan)  (t.gl  IKg/MHp)  (»i/i») 
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Figure  9.  Sample  Test  Report 
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Figure  13.  019  2"  Band  Width  Comparison 


46 


C 
(0 
CD 

CM 

•C 

o 

•** 

m 

CM 

CO 

o 
o 




O" ™Q" Q"" 

* 

o 

O     Q      o       o 

0 

0 

o 

•        * 

U-   UL 

<  m 

C\J   C\J 

a  ° 

CO  CQ 

•  o 

C\J 


CD 


O) 

c 
o 

CO 

0) 

'55 
o 
a. 

E 
o 

O 


CO 


,d  -  ejnuBj  jo  AuiiqeqoJd  linqieM 


sz 
■o 

C 
(0 

m 


u 

CQ 


CM 

CO 

o 




O  ©      ©     G 

O     O     O      O       ( 

>        6           A 

A 

<  5 

J=  -fi 

CQ   CQ 

•  o 

a) 


en 

c 
a> 


ca 

a>        "^ 

'55 
o 
a. 

E 
o 
o 


cvj 


CD 


,d  -  ejniiej  io  Auiiqeqojd  linqie/w 

Figure  14.  008  2"  Band  Width  Comparison 
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Figure  15.  008  10"  Band  Width  Comparison 
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Figure  16.  019  2"  Weibull  Model-Merged  Composite  Data  Comparison 
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Figure  17.  019  10"  Weibull  Model-Composite  Data  Comparison 
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Figure  18.  008  2"  Weibull  Model-Merged  Composite  Data  Comparison 
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Figure  19.  008  10"  Set  A  Weibull  Model-Merged  Composite  Data 

Comparison 
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Figure  20.  008  10"  Set  B  Weibull  Model-Merged  Composite  Data 

Comparison 
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Figure  21.  008  10"  Composite  Band  Width  Comparison 
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Figure  22.  008  (All)  10"  Weibull  Model-Merged  Composite  Data 

Comparison 
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Figure  23.  019  Composite  Data-Metric  Relationship 
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Figure  24.  008  Merged  Composite  Data-Metric  Relationship 
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Figure  25.  019  Composite  Data  (per  fiber)-Metric  Relationship 
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Figure  26.  008  Composite  Data  (per  fiber)-Metric  Relationship 
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Figure  27.  008-019  Composite  Strength  Relationship 
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Figure  28.  019  Fiber  Data  Band  Width  Comparison 
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Figure  29.  008  Fiber  Data  Band  Width  Comparison 


62 


G> 

3 

'co 

u. 

»*- 
o 

>» 


n 
o 


o 

a 

iZ 
© 


E 


CO 

c 
o 


,d  •  ejniiBj  *o  Ajinqeqojd  linqieM 


Figure  30.  019  Weibull  Model-Merged  Fiber  Data  Comparison 
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Figure  3 1 .  008  Weibull  Model-Merged  Fiber  Data  Comparison 
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Figure  32.   019  Fiber  Data-Metric  Relationship 
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Figure  33.  008  Fiber  Data-Metric  Relationship 
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Figure  34.  008-019  Fiber  Strength  Relationship 


67 


019   Fiber/Composite   Probability  of  Failure 


Li- 

I 

<D 

i— 

'(5 
Li- 
's 


CO 

n 
o 


3 


100 


Strength  (per  fiber)  -  [gm] 


Figure  35.  019  Fiber/Composite  Probability  of  Failure 
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Figure  36.  008  Fiber/Composite  Probability  of  Failure 
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Figure  37.   008  Composite  Reliability  Prediction 
Based  on  Fiber  Statistics 
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APPENDIX   B 


TABLE  I.       CALIBRATION  DATA 


Calibration 
Load  (lbs) 


Test  Instrument 
Load  (lbs) 


Calibration 
Drift  (%) 


Before 


During 


After 


0 

0.059 

0.059 

0.059 

— 

11.023 

11.067 

11.008 

11.127 

+.54 

22.046 

22.076 

22.017 

22.135 

+.27 

33.069 

33.084 

33.025 

33.202 

+.36 

44.092 

44.092 

43.974 

44.092 

+0.0 

55.116 

55.101 

55.041 

55.160 

+.11 

66.139 

66.109 

66.109 

66.168 

+.09 

77.162 

77.117 

77.117 

77.236 

+.15 

88.185 

88.185 

88.125 

88.244 

+.07 

99.208 

99.193 

99.134 

99.252 

+.06 

110.231 

110.200 

110.140 

110.200 

+0.0 
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TABLE  II.      AS4-008  2  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 

FAILURE  LOAD  EXPERIMENTAL1  STRAIN  AT 

(Kg)  CERTAINTY  FAILURE 

(mm) 

50.44  1  .0330 

58.52  0  .0374 
56.54  1  .0405 
51.76  0  .0410 
58.79  0  .0362 
55.09  0  .0355 
53.99  0  .0350 
53.88  0  .0351 
57.66  0  .0373 
55.84  0  .0357 
56.48  0  .0385 
58.71  1  .0390 
57.61  0  .0362 
55.25  0  .0359 
55.30  1  .0361 
58.25  0  .0388 
55.17  0  .0358 

57.53  0  .0375 
54.63  0  .0358 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  II.      AS4-008  2  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 
(CONT) 

FAILURE  LOAD                                      EXPERIMENTAL1  STRAIN  AT 

(Kg)                                                      CERTAINTY  FAILURE 

56.00  0  .0355 

55.01  0  .0352 
53.66  1  .0425 
54.25  0  .0452 
55.28                                                            0  .0367 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  III.     AS4-008  10  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA  (SET  B) 


FAILURE  LOAD  EXPERIMENTAL1  STRAIN  AT 

(Kg)  CERTAINTY  FAILURE 

(mm) 


55.89  0  .1710 

53.99  0  .1799 

53.77  1  .1732 

55.68  0  .1856 

55.73  0  .1759 

54.07  1  .1748 

54.98  1  .1779 

53.10  1  .1702 

55.33  0  .1782 

54.58  0  .1770 

47.49  1  .1602 
56.24  0  .1773 
55.60  1  .1788 
51.57  1  .1658 
53.83  1  .1761 
54.17  0  .1776 
54.66  1  .1773 
35.57  0  .1181 
51.87  0  .1645 
52.81  0  .1710 
55.68  1  .1809 
54.44  1  .1753 
56.62  0  .1794 
54.17  1  .1714 

53.50  0  .1717 
53.96  1  .1733 
53.05  1  .1735 
52.56  0  .1688 
54.31  1  .1745 
53.15  1  .1763 

Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  IV.      AS4-008  16  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 


FAILURE  LOAD  EXPERIMENTAL  STRAIN  AT  FAILURE 

(KG)  CERTAINTY1  (mm) 


53.42  0  .2745 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  V.      AS4-008  5  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 


FAILURE  LOAD  EXPERIMENTAL  l  STRAIN  AT  FAILURE 

(Kg)  CERTAINTY  (mm) 

55.52  0  .0896 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  VI. 

AS4-019  2  INCH  C 

AND 

DISPLACEM 

AJRE  LOAD 

EXPERIMENTAL 

(KG) 

CERTAINTY1 

54.50 

0 

57.74 

0 

58.79 

0 

60.24 

0 

58.93 

0 

54.98 

0 

51.60 

0 

62.95 

1 

47.25 

0 

52.08 

0 

52.30 

0 

53.29 

0 

51.57 

0 

53.74 

0 

55.03 

0 

53.91 

1 

51.89 

0 

54.58 

0 

54.79 

0 

51.03 

0 

55.41 

0 

STRAIN  AT 

FAILURE 

(mm) 

.0365 
.0339 
.0328 
.0333 
.0372 
.0336 
.0341 
.0338 
.0320 
.0344 
.0362 
.0361 
.0361 
.0345 
.0355 
.0347 
.0337 
.0403 
.0384 
.0332 
.0447 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  VI.      AS4-019  2  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 
(CONT) 

FAILURE  LOAD  EXPERIMENTAL  STRAIN  AT 

(KG)  CERTAINTY1  FAILURE 

(mm) 

51.95  0  .0397 

54.76  1  .0448 

50.09  1  .0330 

51.46  0  .0349 

50.58  1  .0339 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  VII.        AS4-019  10  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA 

INTRINSIC  LOAD  EXPERIMENTAL  l 

(Kg)  CERTAINTY 

44.89  0 

42.38  0 

44.86  0 

47.15  0 

42.18  0 

50.60  0 
51.91  0 
49.08  0 
48.31  0 

43.47  0 
47.69  0 

45.61  0 

46.79  0 
50.35  0 
47.29  0 
48.44  0 
47.54  0 
45.18  0 

50.80  0 

44.48  0 
52.22  0 
46.53  0 
40.80  0 
48.05  0 

Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  VIII.       AS4-008  10  INCH  GAUGE  LENGTH  LOAD 
AND  DISPLACEMENT  FAILURE  DATA  (SET  A) 


INTRINSIC  LOAD 
(Kg) 

53.19 

45.85 

57.52 

43.24 

49.40 

53.66 

54.03 

54.28 

43.71 

53.51 

51.73 

50.67 

48.65 

48.62 

54.70 

50.00 

51.34 

51.05 

53.32 

51.47 

45.82 

44.04 

45.63 

48.25 

48.74 


EXPERIMENTAL  1 
CERTAINTY 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Note  1 :  '0'  -  No  known  experimental  artifacts 

T  -  Known  experimental  artifact  -  intrinsic  strength  >  failure  strength 
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TABLE  IX.     WEIBULL  DISTRIBUTION  PARAMETERS 
Test  Sample         a  p         AL     TLM    AM      BM     TMU    AU1 


008Fiber(lcm) 

— 

-- 

019Fiber(lcm) 

— 

— 

008  Strand 

33.62 

57.07 

(2") 

008  Strand 

15.22 

51.81 

10"  (Set  A) 

008  Strand 

26.00 

56.09 

10"  (Set  B) 

008  Strand 

15.83 

54.53 

(Merged) 

019  Strand 

14,75 

56.52 

(2") 

019  Strand 

17.48 

48.34 

(10") 

5.23      3.75       5.23     24.76    32.00      5.23 
3.31       9.62       5.08     23.75     30.00      5.08 


Note  1:  Strand  distribution  is  presented  in  a  unimodal  format. 
Fiber  distribution  is  presented  in  a  trimodal  format  to 
facilitate  data  entry  into  the  LLS  program.  Abbreviations  are 

A  -  a  B  -  P  T  -  transition  point 

L  -  Lower  M  -  Middle  U  -  Upper 
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APPENDIX  C 

I.     COMPLIANCE  CALCULATIONS 

To  obtain  the  most  accurate  load  versus  displacement  representation 
possible,  the  raw  cross  head  displacement  data  must  be  adjusted  to  take  into 
account  the  compliance  of  the  mechanical  system  itself.  The  displacement 
recorded  (3)  during  data  collection  is  a  combination  of  the  compliance  of 
the  mechanical  system  (d')  and  that  of  the  test  strand  (9")  such  that 

d  =  d'  +  d"  (1) 


Knowing 


9,  =  P/.  (2) 


'k 


where  P  is  the  load  and  k  '  is  the  equivalent  spring  constant  of  the 
mechanical  system  (which  is  a  constant).  For  the  strand  itself, 

a  =  Ee    or    P/=E —    or  9"  =  —  (3) 

/A  /  AE 

where  A  is  the  cross-sectional  area  of  the  strand 
E  is  Youngs  Modulus 
/  is  the  gauge  length  of  the  sample 

Note  that  since  the  test  strands  are  considered  identical,  1/AE  is  a 
constant. 
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By  substituting  equations  (2)  and  (3)  into  (1)  for  two  samples  of  different 
gauge  lengths,  the  resulting  system  spring  constant  is 


k'  = 


h 

->1    ' 

9,/2 

aaz, 

_  Pi 

p2_ 

(4) 


Since  in  data  interpretation  it's  desirable  to  express  k'  as  a  function  of  a 
load  (P),  we  take  the  case  where  the  load  is  the  same  for  each  data  point, 
thus  Pi  =  P2  =  P,  therefore  equation  (4)  simplifies  to 


k  '  = 


P(/r/i) 


(5) 


Knowing 


A3'  = 


AP 


*'(P) 


we  can  solve  for  d'  at  any  point  by  integrating,  thus  resulting  in 


a\  = 


/ 


*'(P) 


dP 


By  using  equation  (5),  the  system  spring  constant  (k  ')  and  compliance 
(8'),  both  functions  of  load,  were  calculated  at  each  point  that  a  common 
value  of  the  load  occurred  during  the  testing  of  two  separate  gauge  length 
samples.  Ideally,  compliance  data  should  be  taken  with  a  constant  load 
control  system  instead  of  a  constant  cross-head  rate  control  system  for 
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direct  data  reduction.  This  feedback  control  configuration  is  not  available 
for  the  testing  system  used  in  this  investigation;  as  an  alternative,  the  data 
from  two  different  gauge  length  samples  were  manually  reduced  and 
common  loads  identified  in  order  to  conduct  the  analysis.  Another 
approach  that  could  be  taken  would  be  to  curve  fit  the  test  data,  establish  an 
equation  (displacement  as  a  function  of  load)  for  each  set  of  data,  and  then 
simply  calculate  the  displacements  at  the  desired  loads. 

Four  samples  of  varying  gauge  lengths  (2,  5,  10,  and  16  inches)  were 
tested  and  analyzed  to  obtain  the  system  compliance.  The  samples  were 
produced  as  described  in  the  Background  section  using  AS4-008  composite 
strands.   All  combinations  of  d'(P)  computed  from  the  four  samples  were 
compared  and  the  resulting  graphical  comparison  is  presented  in  Figure  1 . 
As  shown,  system  compliance  closely  correlates  in  three  (2-5  inch,  2-10 
inch,  and  2-16  inch)  of  the  six  combinations.  Since  the  calculated  value  of 
D'(P)  at  low  loads  for  the  2-5  inch  gauge  length  compliance  is  negative,  the 
final  compliance  calculation  excludes  this  result.  Thus,  the  final  system 
compliance  was  derived  by  averaging  the  2  and  16  inch  gauge  length 
compliance  data  with  the  2  and  10  inch  gauge  length  compliance  data  which 
is  presented  in  Figure  2.  The  equation  for  the  system  compliance  was 
obtained  by  using  a  third  order  polynomial  curve  fit  through  the  averaged 
data  and  is 

3'(P)  =  9.3062E-5  +  2.8145E-5  P  +  7.9643E-8  P2  +  4.7468E-10  P3 

It  should  be  noted  that  this  equation  is  valid  only  for  loads  up  to   145 
lbs  since  extrapolation  of  the  curve  fit  polynomials  are  not 
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necessarily  valid  beyond  this  point.      However,  this  is  sufficient  since 
the   maximum   loads   observed   during   testing   do  not  exceed  this   limit. 
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Figure  1.    System  Compliance  Comparison 
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APPENDIX  D 

I.    MERGING  OF  STRENGTH  DATA  FROM  DIFFERENT 

SAMPLE  DIMENSIONS 

There  are  two  different  mathematical  normalizations  of  the  strength 
data  measured  from  samples  of  different  dimensions  (in  the  context  of  this 
investigation,  different  gauge  lengths).  The  underlying  assumption  of  the 
mathematical  normalization  is  that  the  failure  is  a  extreme  value  process 
(weakest  link)  and  that  the  strength  of  each  segment  of  the  sample  (the  link) 
is  Independent  and  Identically  distributed. 

When  data  are  represented  in  a  transformed  Weibull  cumulative 
probability  of  failure,  a  horizontal  shift  of  the  data  can  be  made  which,  in 
effect,  shifts  the  location  parameter  of  the  assumed  two  parameter  Weibull 
distribution  model.  This  results  in  a  shift  of  the  lower  tail  of  the  Weibull 
probability  curve  down.   The  derivation  for  this  correction  of  gauge  length 
is  presented  in  Appendix  C  of  Ref.  7.  It  should  be  noted  that  this  method 
can  only  be  utilized  on  a  two  parameter  Weibull  distributed  function 
(which  is  linear  in  the  Weibull  probability  of  failure  coordinates)  or  a 
three  parameter  function  that  can  be  linearized  into  a  two  parameter 
function. 

A  vertical  shift  of  the  data  can  also  be  made  which  maintains  the  same 
location  parameter  but  shifts  the  Weibull  probability  curve  up  or  down  on 
its  vertical  scale  in  order  to  account  for  the  change  in  probability  associated 
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with  a  change  in  the  metric  dimension.  It  can  be  shown  that  the  reliability 
of  a  chain  (fiber)  is  equal  to  the  reliability  of  a  link  in  the  chain  (the  metric 
dimension)  raised  to  the  power  of  the  number  of  links  in  the  chain  (n),  or 
in  terms  of  probability  of  failures 

(1  -  Fn(L))  =  (1  -  F(L))n 
where     Fn(L)   is  the  probability  of  failure  of  the  chain  and 

F(L)      is  the  probability  of  failure  of  a  link  in  the  chain 

This  assumes  that  the  links  of  the  chain  are  independent  and  identically 
distributed.  The  equation  can  be  written  as 

ln(l  -  Fn(L))  =  n  In  (1  -  F(L))  or 

ln(-ln(l  -  Fn(L)))  =  In  n  +  ln(-ln  (1  -  F(L))) 

which  is  observed  to  be  in  the  form  of  x'  =  c  +  x.  Therefore,  the 
distribution  of  the  chain  preserves  any  and  all  distributions  of  the  link  and 
the  distribution  of  the  chain  is  simply  a  vertical  shift  of  the  link  by  In  n  (the 
physical  size  of  the  chain).  Note  that  the  use  of  ln(-ln(l  -  Fn(L)))  implies  a 
series  of  weakest  links  (Type  HI,  extreme  value  statistics)  and  if  it  is  linear, 
then  it  implies  a  Weibull  distribution.  The  same  procedure  used  in  the 
above  mathematical  derivation  is  also  valid  for  a  composite  with  load 
sharing.    Applying  these  procedures  results  in 

ln(-ln(l  -  Hmn(L)))  =  ln(mn)  +  ln(-ln(l  -  W(L))) 
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where  Hmn(L)  is  the  probability  of  failure  of  a  composite  of  length  / 

(given  mn)  where  /  is  the  gauge  length 
W(L)     is     the  probability  of  failure  of  a  composite  at  the 
desired  metric  length 
n        is    the  number  of  fibers  in  the  bundle 
m        is    the  number  of  bundles  (segments  of  the  length 
of  the  metric)  in  the  gauge  length 

Note  that  this  again  is  in  the  form  of  x'  =  c  +  x  and  that  the  transformation 

simply  involves  a  shift  in  the  vertical  axis.    Rearranging,  we  have 

ln(-ln(l  -  W(L»)  =  ln(-ln(l  -  Hmn(L)))  -  ln(mn) 

This  relation  is  used  to  transform  strength  data  measured  from  different 
gauge  lengths  into  one  common  gauge  length.  This  normalization  of  all 
data  into  one  reference  metric  dimension  allows  the  merging  or  pooling  of 
the  aforementioned  strength  data  for  comparison  with  the  analytical  Local 
Load  Sharing  probability  predictions. 
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APPENDIX   E 
I.    MAXIMUM  LIKELIHOOD  ESTIMATOR 

A. BACKGROUND 

The  Maximum  Likelihood  Estimator  (MLE)  method  is  a  method  of 
estimating  the  parameters  of  a  model  (curve)  for  a  given  set  of  data.  There 
are  several  advantages  of  the  MLE  method  over  the  least  squares  method. 
The  least  squares  method  assumes  equal  weighting  of  the  data  in  the 
linearized  domain  and  does  not  account  for  data  clustering.  In  addition,  it 
can  only  be  used  with  exact  (x  =  xj  )  data,  where  x  is  the  true  realization  of 
the  data  and  x[  is  the  experimental  data.  The  MLE  method,  on  the  other 

hand,  correctly  weights  the  data  by  the  probability  of  occurrence  and  it  is 
applicable  for  exact  (x  =  xj  ),  interval  (xj  <  x  <  xj+i),  or  censored 
(x  >  xj )  data.    A  data  set  may  consist  of  one  or  more  of  the  three  classes 
of  data.  The  MLE  method  is  statistically  sufficient,  asymptotically 
consistent,  asymptotically  efficient ,  and  converges  to  the  underlying 
parameters. 

The  principle  behind  the  MLE  method  is  to  determine  the  model 
parameters  of  the  assumed  model  that  has  the  maximum  likelihood,  or 
probability,  of  repeating  a  particular  set  of  data  from  all  possible 
repetitions.  In  developing  the  mathematical  representation  of  the  MLE 
method,  it  is  assumed  that  the  combined  probability  of  all  the  data 
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occurring  during  testing  is  the  product  of  the  probability  of  each  data  point 
in  the  set,  or 


n    f(xi,0)=  f(Xl,0)   f(x2,0)   f(x3  ,0) 
i=  1 


f(xn,0) 


(1) 


where  f(xi50)  is  the  probability  of  one  of  n  data  points  given  the  model 

parameters  0.  This  assumes  that  the  data  were  collected  from  identical, 
independent  samples.  The  likelihood  (L)  is  therefore  defined  as  this 
combined  probability,  or 


L=    n    f(xif9) 
i=l 


(2) 


One  notes  that  to  solve  for  the  likelihood  (L),  the  probability  function 
fCxj.G)  must  be  known,  which  includes  both  the  data  and  the  model 

parameters.  The  problem,  of  course,  is  that  the  parameters  are  not  known. 

One  approach  to  solving  this  problem  is  to  conduct  an  iterative 
calculation  process  such  that  you  assume  values  of  a  and  p,  and  solve  for 
(L).  The  likelihood  (L)  is  then  plotted  as  a  function  of  a  and  p,  resulting 
in  a  3-D  plot.   Different  values  for  a  and  p  are  then  chosen  and  another  L 
is  calculated  and  plotted.  This  process  is  continued  until  the  plot  can  be 
analyzed  and  the  peak  or  maximum  likelihood  determine.  These  values  of 
a  and  p  become  the  estimated  parameters  of  the  model. 

There  are  several  problems  associated  with  this  process.  First,  since 
f(xi50)  <  1  ,  it  implies  that 
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n 

n  f(xi,9)«i 

i  =  1 
and  for  large  sample  sizes,  the  computer's  numerical  capacity  can  be 
underflowed.   In  addition,  the  evaluation  of  where  the  maximum  likelihood 
occurs  involves  the  use  of  a  concept  dealing  with  Confidence  Regions, 
which  turn  out  to  be  coupled  with  associated  Confidence  Intervals  for  the 
parameters.  The  evaluation  of  (L)  can  then  become  quite  complicated. 

Another  solution  to  the  problem  is  to  solve  directly  for  a  and  p.  This 
solution  maximizes  the  likelihood  (L)  analytically  by  taking  the  derivative 
(with  respect  to  the  model  parameters)  of  equation  (2)  and  setting  them 
equal  to  zero.   For  the  two  parameter  Weibull  distribution,  these 
mathematical  steps  can  be  performed  explicitly  resulting  in  two  equations 
with  the  two  unknown  model  parameters  (a  and  p)t  These  two  equations 

are  then  solved  simultaneously  for  the  model  parameters.  As  it  rums  out, 
the  equations  for  a  and  p  are  decoupled  such  that  one  parameter  (normally 
a)  can  be  optimized  and  then  the  other  parameter  (p)  calculated. 

B.  CENSORED  DATA  WITH  KNOWN  BUT  DIFFERENT 
LOWER  BOUNDS  (WEIBULL  DISTRIBUTION) 
1.     Definition   of  Variables 

n   =  Total  number  of  tests 

m  =  Number  of  tests  of  known  realized  strength  (no  uncertainty) 
Xj  =  Exact  data 

xj  =  Realized  data  where  \{<  x;  (uncertain  data) 
0    =  Parameters  of  the  assumed  model 
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The  likelihood  is  the  product  of  all  the  individual  probabilities,  but  the 
exact  value  for  the  uncertain  data  is  unknown.  What  is  known  is  that  the 
exact  data  is  at  least  as  great  as  or  equal  to  the  uncertain  data.  Therefore 
the  likelihood  is 


L  = 


~ 

m 

n  f(Xi,e) 

\f7"  f(x,6)dx 

f*T  f(x,0)dx 

X2 

i=  1 

L      ' 

or 


L  = 


m 


n  f(xbe) 

i=  1 


n 


i  =  m 


X  4x'e) 


dx 


Knowing  that 


f(x,0)dx  =  1  -   (  X'f(x,6)dx  =  1  -  F(Xi,9) 

Jo 


the  likelihood  is  then 


L  = 


m 


n  f(xite) 

i=  1 


n  (i-F(Xi,e)) 


i  =  m 


For  a  2  parameter  Weibull  distribution  we  know 


f(x,9)  = 


ax 


P 


a-l 


exps 


Ipi 


and 
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1-F(xi,0)  =  exp\-P] 


Therefore 


or 


L  = 


m  a-l 

„       axi  i     /xi 

n     exp\  - 


a 


L  = 


i=l       p 


m 

m    -ma 

a  p       n 

i=  1 


\\T 


n    exp\-  — 


Xi 


(a-l) 


exp 


a 


n  r^ 

x 
n     exp 


1  =  m 


IP/ 


2.     Maximizing  L  with  Respect  to  a  and  p 

To  remove  the  multiplicative  products  from  the  above  equation, 
we  take  the  natural  log  (In)  of  L  such  that  ln(L)  =  L.   Therefore 


-a 


•a     "      ^a 


L  =  m  ln(a)  -  ma  ln(p)  +  £  (a  -  1)  ln(xj)  -  p     £  x"  +  P   '  £  xi 

i  =  1  i  =  1 


i  =  m 


To  maximize  L  (equivalent  to  maximizing  L)  with  respect  to  a  and  p  ,  the 
partial  derivatives  are  taken  and  set  equal  to  zero,  or 


^  =  0   and   ^  =  0 

5a  5p 


This  results  in 

5L  _1      x^ 

—  =  0  =  ma     +  X  ln(xi) "  ln(P) 

5a  i  =  l 


m-p 


-a 


m  n        ry 

Lx'  +  Lxi 

i  =  1         i  =  m 


-|\ 


-a 


-1/ 


^(xTlnCxj)) 


i  =  1 
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n 

~<x       ^ 

I 

Xj    ln(xj) 

i  =  l 

and 


—  =  0  =  -ma  +  ap     Zxi+2>i 
5(3  vw        ™ 


The  above  equation  simplifies  to 


(3) 


m=p 


-a 


m  n         a 

2,x'  +  2.xi 

i  =  1  i  =  m 


or 


P  =  /m- 


m  n        n 

2,*'  +    ZXi 
i  =  1  i  =  m 


(4) 


Substituting  (4)  into  (3)  and  solving  for  a  and  p  produces 


i    m  /  \        n  l^a    *~ 

—     y,(XiIln(Xi)J  +  —  S\Xi    lll(Xi 


i  =  1 


i  =  1 


a  = 


1      ">,       a        1      "~a 

—  L«i  +  —  L«i 

m  i=i  m  i=m 


-Iln(Xi) 
m  i=i 


and 


m  i       n  a 

K^X^  +  mX** 

1  =  1  i  =  m 


a 
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We  note  that  the  shape  and  scale  parameters  a  and  p  are 
decoupled  but  a  is  implicit.  Either  an  iteration  calculation  or  numerical 
methods  are  used  to  solve  for  the  shape  parameter  a.  Once  a  is  found,  (3  is 
calculated  directly.    We  further  note  that  the  equations  for  a  and  (3  can 
also  be  used  for  non-censored  data  by  simply  letting  m  =  n. 
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